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, We review the state-of-the-art of our understanding of heavy quark fragmentation. Recent 

e + e~ data for B mesons are compared to the most up-to-date theoretical predictions, and 
the need for inclusion of a non-perturbative component is discussed. Experimental analyses 
in moments space are suggested, and it is pointed out how perturbative and non-perturbative 
5^ ' contributions are to be properly matched. Failure to do so can result in large phenomenological 

discrepancies. An example is given for B + hadroproduction at the Tevatron. 



1 Introduction 




This talk will be devoted to a review of our present understanding of heavy quark fragmentation, 
i.e. of the processes where a heavy quark is produced in a hard collision and then hadronizes 
before decaying. For the quark to be "heavy", its mass has to be larger than the QCD scale 
Aqcd — 200 — 300 MeV. The top quark, however, is too heavy and decays weakly before 
hadronizing. Throughout our discussion "heavy quark" will therefore have to be understood as 
a charm or a bottom quark. 

Let us start by considering, for the sake of definiteness, an experiment observing a B meson. 
Since we know it to contain a b quark, whatever the initial state of the process we can envision 
its production to proceed as roughly sketched in Fig. [I]: A so-called hard process, described by 
perturbative QCD (pQCD), produces a bottom quark. Subsequently, soft strong interactions 
(the non-perturbative 'np' blob) bind this heavy quark to light ones into a hadronic state, 
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Arbitrary separation 




Figure 1: Schematic view of how the production process of an observed B meson can be separated into pertur- 

bative and non-perturbative contributions. 

resulting into the observed B meson with only a fraction z of the original momentum of the b 
quark. 

From such a simplified description alone, it appears clear that any separation between "hard" 
and "soft" processes, which is to say between perturbative and non-perturbative QCD, is at best 
arbitrary. We shall focus more on this issue in the following. 

People knowledgeable with light quark fragmentation only might not be immediately famil- 
iar with the following fact: In the heavy quark case the pQCD term does not need collinear 
factorization to be performed. It is the large mass m of the heavy quark that takes care of 
regulating the collinear singularities, acting as an infrared cutoff. Hence, pQCD alone can yield 
a finite result and be predictive at the 'leading twist' level. Further discrepancies with the ex- 
perimental results (aside, of course, unknown higher order perturbative contributions) can then 
be attributed to a non-perturbative part, to be studied in the form of power corrections. 

Many analyses, and also intuitive arguments, predict for the leading term of such power 
corrections a A/m behaviour, A being a hadronic scale of the order of a few hundred MeV. These 
corrections are therefore parametrically small, m being larger than A, but can be numerically 
large. For instance, they are definitely larger than typical power corrections in event shapes 
studies at LEP, where the suppressing scale is ~ 90 GeV rather than the heavy quark mass. 
Heavy quark fragmentation is therefore an ideal place to observe and analyze such corrections. 

2 Perturbative Issues 

Let us consider the production of the heavy quark Q in e + e~ collisions, e + e~ — > j/Z —* Q + X. 
An experimentally observable variable usually looked at is the energy fraction (with respect to 
the beam energy) of a heavy hadron H containing the heavy quark, xe = Eh / E^ earn . 
Fixed order perturbative calculations 

Hi 

are available up to order a 2 s . These results do 
however contain large logarithmic terms which need to be resummed to all orders: 

• log(s/m 2 ) terms are large when the centre-of-mass energy yfs is much larger than the 
heavy quark mass m. Such a situation is easily met at LEP energies for both charm and 
bottom production; 

• 1/(1 — xe) and log(l— xe)/{1 — xe) terms are large when the energy of the observed particle 
is close to the maximum allowed one, and are due to gluon radiation being inhibited close 
to the phase space boundaries. 

Mele and Nason@ first considered the resummation of log(s/m 2 ) terms up to next-to-leading 
logarithmic (NLL) level. They achieved it be rewriting the e + e~ "f/Z — > Q + X differential 



cross section in a factorized form]], 



a N (V~s, m) = C N {^s, ii F )B$(ji F ,m) + 0{{m/^sf) . (1) 

The factorization scale fj, F separates the two functions C and D im , and Altarelli-Parisi evolution 
can be used to resum the collinear logarithms in D mi (u, F ,m) by evolving from an initial scale 
Hqf ~ m up to the hard scale \x F — \/s, so that £)™(^,m) = En(/j, f , fj,Q F ) D 1 ^ 1 (fiQ F , m) . 

It is to be noted that all the information concerning the hard scattering process is now 
contained in the coefficient function Cn{V~s, /j, f ). The initial condition D 1 ™ (u. 0F , m) , on the 
other hand, only contains information about physics taking place around the heavy quark mass 
scale m. This process independence of the initial condition function has been recently established 
in a more formal setting!. It is worth reminding that, while process-independent, D™{no F ,m) 
is of course still factorization scheme dependent, like the coefficient function and the Altarelli- 
Parisi evolution factor. Both Ref J and Ref.! make use of the MS scheme. 

After performing the resummation of the log(s/m 2 ) terms one still has to take care of the 
potentially terms related to suppression of radiation close to the phase space borders. Such 
terms show up as large logarithms in the large- N limit of the Mellin transforms. For instance, 
in this limit the initial condition reads: 
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Resummation for these so-called Sudakov logarithms was performed at the leading log (LL) 
level! and at NLL level, but in a process dependent way!. In Ref.!, on the other hand, NLL 
resummation has been revisited in a fully process independent way, exploiting the factoriza- 
tion properties of D 1 ^ 1 (^iQ F ,m) and providing NLL resummed expressions for both the e + e~ 
coefficient function and the process independent initial condition. 

The resummation of all the large logarithms present in the perturbative calculation for 
DN(y/si m ) allows for a prediction which is at the same time reliable and accurate. The depen- 
dence on the unphysical factorization scales fj, F and [jlq F , for instance, is greatly reduced and 
under control after the Sudakov logarithms are resummed with NLL accuracy ti. 

3 Non-Perturbative Issues 

Despite having achieved a reliable perturbative prediction, inclusion of non-perturbative effects 
is still mandatory for a meaningful comparison with the experimental data, for instance the 
ones recently published by the ALEPH0 and SLdH Collaborations. After calculating a cross 
section for, say, b quark production, the one for B meson is obtained by convoluting it with a 
non-perturbative component: 

D B N = D b N D% . (3) 

Pretty much at odds with the available literature, I shall consider comparisons between 
theory, and experiments in iV-moments space rather than with differential distributions in x F 
space @M 

The most obvious reason for doing so is that only in moments space a smooth departure 
of the experimental results from the perturbative prediction can be seen: low-TV moments are 



b From now on we shall make use of Mellin moments, Djv = x N 1 D(x) d(x), which turn convolutions into 
products. 




closest, and the gap widens as N increases. Indeed, the non-perturbative contribution can be 
predicted to behave like 00 

0% = l-(N-l)± + o(^) . (4) 

In xe space, on the other hand, one usually tries to compare a whole curve, and the shape of 
the perturbative prediction has clearly little to do with the experimental one, especially in the 
xe — 1 region, where the peak is. This is due to the fact that all-order power corrections become 
important here, and prompt one to include from the very start a "large" non-perturbative term, 
for instance in the form of some smearing function like the Peterson et al. one, in order to fit 
the data. Doing so, however, immediately obscures the fact that non-perturbative contributions 
can be seen as a small A/m correction to the perturbative result. 

Figure |3| clearly shows this point. The purely perturbative result fails in describing even the 
lowest N moment (see left panel), hence showing that inclusion of a non-perturbative component 
is mandatory. However, its job looks much worse in xe space (right panel), particularly in the 
region around and beyond the peak. But it suffices to fit a 1-parameter non-perturbative form| 
to the Z?2 = (xe) point in N space, to produce a curve which describes very well the low-iV 
region, and even xe data much better (albeit not perfectly). We wish to stress that no effort has 
been made in this case to achieve a particularly good description of the whole xe distribution, 
the emphasis being rather on showing that a "reasonable" functional form, fitted to a single point 
in moments space in terms of a numerically small A/m ~ 0.1 correction, can already produce a 
decent agreement. Good fits of the low-iV moments, on the other hand, will be important for 
the issue we are going to discuss in the next Section. 



4 Phenomenology 

The theoretical machinery of fragmentation functions for heavy quark can be used to make 
predictions for processes other than e + e _ production, exploiting the property of process inde- 

c In this case a D np (x;a) = (a + l)(a + 2)x a (l — x) form (or, rather, its Mellin transform) was used, and the 
fit returned a = 27.45. It can easily be seen that this form is compatible with the leading power correction (^) 
upon replacing a with 2m/ 'A. This equality suggests A ~ 350 MeV, in line with our intuitive expectations. 
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Figure 3: On the left, the experimental data from CDFliJ on B + production in pp collisions at the Tevatron, 
compared to theoretical predictions. On the right, moments of ALEPH u data on B fragmentation in e + e~ 
collisions compared to a purely perturbativerfcJpeoretical prediction and to various models/fits for non-perturbative 
contributions. Plots from Ref. EJ. No Sudakov resummation was included in this case. 



pendence of the initial condition D mi . For instance, bottom cross sections at large transverse mo- 
mentum rtr in hadronic collisions can be calculated, providing a resummation of large \og{pT/m) 
terms EIH At the same time, inclusion of the non-perturbative component determined from 
e + e~ fits allows to make predictions for the production of B mesons, which are the particles 
which can be directly observed. 

Obtaining accurate theoretical predictions requires however great care in at least two in- 
stances: 

• a "non-perturbative term" is not an observable quantity. It cannot be determined in 
absolute terms, but only relatively to how one defines the perturbative part and its pa- 
rameters. Therefore, when fitting, say, e + e~ data one extracts a non-perturbative function 
which should then be used only together with a perturbative description of the same Jcind 
(leading, next-to- leading, resummed, etc.) as the one it has been determined with, and 
with the same parameters {KqcDi m, ■■■ )'■, 

• different processes may depend on different details of the non-perturbative contribution. 
For instance, the peak in xe space is clearly the most prominent feature in e + e~ processes, 

isions, 
. This 



but not the most important one when calculating B meson production in pp col 
where instead a good determination of the moments around N ~ 4 is crucial 11311^ 
is due to the fact that the cross section for production a heavy quark in hadronic colli- 
sions behaves like da b /dpx — A/pj,. Hence, convolution with a non-perturbative function 
returns: 

drr B f A A 

-i—^ dz dp T D n P(z)—5(zp T - PT ) = — D%> , (5) 
dpT J Pr Pt 

with n — 4. 

That these issues are not merely of academic interest is shown by the case of B production 
at the Tevatron 0. Figure |3] shows that the "standard" procedure, used by the experimental 
Collaboration, of evaluating the cross section by convoluting the perturbative calculation for b 
quark with a Peterson et al. function with e = 0.006 underestimates the data by almost a factor 
of three. 

On the other hand, employing a non-perturbative contribution fitted to e + e~ data in mo- 
ments space (the U N = 2 fit"), so as to get a good description of moments around N ~ 4 (and 
taking care of employing the same kind of perturbative description in e + e~! and p«0 processes) 
increases the prediction, bringing it in agreement with the data within the errors 



5 Conclusions 



Heavy quark fragmentation is now a mature subject. Fixed order perturbative calculations 
have been performed up to order a 2 s , collinear and Sudakov resummations are known up to 
next-to-leading logarithmic accuracy. 

Comparisons to experimental data however need the addition of a non-perturbative compo- 
nent. The fairly large contributions expected in this process make it a good one to study power 
corrections. Performing the fits and the comparisons in moments space helps disentangling the 
leading power correction from the higher order ones. 

Factorization of terms related to the hard scattering from terms related to the heavy quark 
mass scale suggests some form of universality for such non-perturbative terms, which can there- 
fore be fitted in one process and used to give predictions in another. 

When transporting this information care must however be taken to properly match the 
perturbative and non-perturbative terms, which are not independently measurable quantities, 
and to carefully consider what details of the non-perturbative function need to be well known. 
Failure to do so may result in predictions with a far larger degree of uncertainty than one might 
expect. 
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